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Abstract 

We point out that, in a class of gauge mediation models using metastable supersymmetry break- 
ing vacua, the minimum of the supersymmetry breaking field in the early universe is dynamically 
deviated from the one in the low energy. The deviation induces coherent oscillations of the super- 
symmetry breaking field, which decays into the gravitinos. For certain parameters, it can produce 
a right amount of the gravitinos to account for the observed dark matter. 



I. INTRODUCTION 

Gauge mediation (GM) of supersymmetry (SUSY) breaking [1] is a natural solution to 
the phenomenological problems such as excessive flavor- changing neutral currents. In spite 
of the successes, building a realistic model is a rather non-trivial task. According to the 
argument of Nelson and Seiberg, an exact U(1)r symmetry is required if the superpotential 
is generic This observation strongly limits possible models, and a lot of efforts have been 
devoted to building a realistic one. 

Recently, Murayama and Nomura has proposed drastically simplified models, focusing on 
metastable vacua in the SUSY breaking sector Even though the entire superpotential 
does not possess the exact U(1)r symmetry, an accidental one exists near a local SUSY 
breaking minimum [4|. Such a scenario is viable as long as the metastable vacua have a 
sufficiently long life time. 

In a class of the metastable SUSY breaking models, there are the local SUSY breaking 
vacua near the origin of the SUSY breaking fields, where the accidental U(1)r symmetry 
exists. On the other hand, the breaking of U(1)r symmetry is necessarily involved in the 
messenger sector to mediate the SUSY breaking to the gauginos in the supersymmetric 
standard model (SSM). Since the SUSY breaking field (denoted by S) linearly couples to 
the messenger fields, the breaking of the U(1)r symmetry induces a linear term of S in the 
Kahler potential. 

Such a linear term forces the SUSY breaking field S to deviate from its minimum in the 
low energy, while the inflaton field dominates the energy of the universe. When the Hubble 
parameter becomes comparable to the mass of S, it starts to oscillate coherently around 
the SUSY breaking vacuum, and then decays into the gravitinos. In the GM models, the 
gravitino is stable and behaves as dark matter (DM) in the universe as long as the R-parity 
is preserved. Thus, DM is generally produced from the SUSY breaking field in a class of the 
GM models using the metastable vacua. In this letter, we study the gravitino production in 
this scheme. 

Before proceeding to the details, let us comment on the previous works. The cosmological 
evolution of the SUSY breaking, field has been studied in the models with the metastable 
SUSY breaking vacua in Refs. [5j. Those literatures assumed so high reheating temperature 
as to make the SUSY breaking sector to be in thermal equilibrium. Then the gravitino is 
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expected to reach thermal equilibrium, and as a result, its abundance exceeds the observed 
DM abundance or it erases the density fluctuation too much, unless the gravitino mass is 
smaller than 16 eV 6|]. For a wide range of the gravitino mass, i.e., from 16 eV to O(10) GeV, 
therefore, the SUSY breaking sector should not be thermalized as long as the standard 
thermal history of the universe is assumed. In this letter, we do not pursue this possibility 
and assume that the SUSY breaking sector never reaches thermal equilibrium. On the other 
hand, Ibe and Kitano discussed the gravitino production from the SUSY breaking field in the 
GM models using the metastable vacua |7j. They assumed a different thermal history taking 
a different set of the parameters; in particular, the reheating temperature in our scenario is 
as high as 10 8-10 GeV, while they considered the low reheating temperature. Furthermore, 
the interaction that shifts the SUSY breaking field from its minimum in the low energy is 
different from that considered in Ref. Q]. 

II. MODEL 

In this section we provide a model of the gauge mediation using the metastable vacua. 
To be explicit, we adopt a model by Murayama and Nomura given in Ref. [3|. The Kahler 
potential and the superpotential for a gauge singlet chiral field S and the messengers / and 
/ are written as 

K = |5| 2 -^ + |/| 2 + |/| 2 , (1) 
W = -ffS + KSff + Mff, (2) 

where the higher order corrections of 0(|5| 6 /A 4 ) are omitted for simplicity in the Kahler 
potential. In the following we take / and / to be in 5 + 5* representation of S77(5)guT) and 
/i 2 , k and M are set to be real and positive without loss of generality. We assign the charges 
U(1)r symmetry as R[S] = 2 and R[f] = R[f] = 0. Then one can see that the messenger 
mass term explicitly violates the U(1)r symmetry to Z2- 

For our purpose we do not need to specify the UV physics above a scale A that provides 
the second term in Eq. ([I]) as well as the first term in Eq. ([2]). We simply note here that there 
are many explicit models that actually lead to this low energy effective theory. (See Ref. [s| 
for examples.) In particular, we do not give the SUSY breaking mechanism explicitly here, 
which is assumed to be such that the first term in Eq. (j2J) is somehow produced. 
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From the Kahler potential and the superpotential given above, one can show that there 
is a SUSY minimum at 

S=--, / = /=-£=■ P) 

K \/K 



On the other hand, SUSY is broken aX S = f = f = 0, which is a metastable local minimum 
as long as 

M 2 > Kfi 2 (4) 

is satisfied, since otherwise one of the messenger scalars becomes tachyonic. Note that the 
second term in Eq. ([T]) produces a positive mass of S, ms = /i 2 /A, around the origin. The 
SUSY breaking scale is dictated by the first term in Eq. (121) , and the F-term of S is given 
by Fs — yU 2 . Requiring a vanishing cosmo logical constant, we can relate the SUSY breaking 
scale to the gravitino mass m 3 / 2 as 



= f V3 m 3/2 M P ) * ~ 2 x 10 9 GeV ( ' , (5) 



where Mp = 2.4 x 10 18 GeV is the reduced Planck scale. The SUSY breaking effects are 
transmitted to the visible sector by the messenger loops. The integration of the messengers 
give rise to the gaugino masses as [j| 

mi ~ 4^1T 1 = 1,2,3. (6) 

Here, 0^1,2,3 and 01,2,3 are the gaugino masses and the gauge coupling constants for 
U(1)y, SU{2)l and SU(3)c in the SSM. We have used the SU(5)gvt normalization for 
the U{1)y gauge coupling constant. We can express k/x 2 /M in terms of the gluino mass m 3 : 

ku 2 / ™~ x _1 

M 

From (J5J) and (171), one can express M as 
We also assume ms ^ A, or equivalently, 

M < A, (9) 

since we consider the dynamics of S (e.g. coherent oscillations and decay), which should be 
described within the low energy effective theory. Using <M), we obtain an upper-bound on 



ioW^vv^-v (7) 

V0.17 VlTeV/ v ; 



m s , 



ms = — ~ 40 lev 



m 3/2 \ / 
lGevJ V 



A UGeV/ V10 14 GeV 



< 2xl0S Gev(i^j-. (10) 

Lastly, let us discuss radiative corrections to the Kahler potential. Integrating the mes- 
senger loop, the relevant corrections are given by 

= k<» + k£\ (ii) 



with 



where we have separated the holomorphic terms and the non-holomorphic ones. Note that 
S is assumed to be much smaller than Mj k so that S sits far away from the SUSY minimum 
(see 03])). As explained in Introduction, the reason why such corrections appear is that the 
i?-symmetry is explicitly broken by the messenger mass term. 

One can check that the radiative corrections (Tl2l) reproduce the result of the Coleman- 



Weinberg potential for S given in Ref. 



3J, up to the order explicitly shown in ( fl2l) : 
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To avoid the mass of S to become tachyonic due to the radiative corrections, we require 

M > ^A, (15) 

throughout this letter. 

As far as the SUSY breaking sector is concerned, the linear term in the Kahler potential 
(fl~3]) does not modify the scalar potential significantly. In the very early universe, however, 
such a linear term makes the minimum of the scalar potential to deviate from the origin. 
Therefore, it is crucial for cosmological evolution of S to take into account the linear term 
in the Kahler potential, as we will show in the next section. 
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III. COSMOLOGY 



Now we consider the cosmological evolution of the SUSY breaking field, S. First let 
us give a sketch how S is deviated from the origin due to the linear term in the Kahler 
potential, and estimate the cosmic abundance. While the F-term of the inflaton dominates 
the universe, the scalar potential of S is approximately given by a 



3# 2 (|S| 2 -^M(S + St) + ... , (16) 



V(S) ~ e K (3H 2 M 2 ), 

5 k 

167T 

where we have assumed that S does not couple to the inflaton in the Kahler potential for 
simplicity b . The scalar potential has a minimum given by 

S c = -|^M. (17) 

167T 

If this minimum S c exceeds A, there is no stable minimum in the low-energy theory in the 
early universe. Depending on the UV theory, the system may settle down at the SUSY 
minimum in this case. To avoid such a situation, we impose S c < A in the following. 

When the Hubble parameter becomes comparable to the mass of S, it starts to oscillate 
around the minimum, S ~ 0, with an initial amplitude S c . The abundance of S is estimated 
as 

n s 3T R ( 5k \ 2 m s M 2 



s 4 V167T 2 / 3m|M|' 

- 1 * " (s) 2 G^r m) doiv) (io4v) ■ w 

where ns is the number density of S, s is the entropy density, Tr denotes the reheating 
temperature, and we have used (jSJ) to eliminate M. We have assumed here that the reheating 
has not completed when H = iris- This assumption is indeed reasonable, since otherwise 
too many gravitinos are produced by thermal scatterings in plasma. 

Several comments are in order. First, we have assumed that the Hubble parameter 
during inflation, Hi, is larger than the mass of S, i.e., Hi > ms- If ms is larger than 



a The linear term of S generically appear in the scalar potential due to the supergrav- 
ity effects 0- One can neglect its effect on the dynamics of S, as long as k > 
0.04(a 3 /0.1) 1 / 2 (lTeV/m 3 ) 1 / 2 (m 3/2 /lGeV) 1 / 2 (A/10 14 GeV). We assume that this inequality is satisfied 
in the following analysis. 

b Even in the presence of the interactions, the following argument does not change qualitatively. 
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Hj, the deviation of S is suppressed by (Hi/rris) 2 . Considering the upper-bound on ms 
given by ffTUl) . however, the assumption Hj > ms is valid except for low-scale inflation 
models. Second, although our model is given only below the scale A, this does not limit the 
application of the above arguments only to the inflation models with Hj < A. For Hj < A, 
the above scalar potential (1161) is obviously valid both during and after inflation. On the 
other hand, for Hi > A, one cannot use (fT6l) during inflation. After inflation, however, the 
Hubble parameter decreases and becomes smaller than A at certain point. If the system 
can be then described by the low energy effective theory, S will quickly settle down at the 
potential minimum ( TT7T) c . 

Next let us consider the decay processes of S. The SUSY breaking field S will decay into 



a pair of the gravitinos, and the decay rate is 



9] 



3/2 " 96vr m 2 M 2 - 32vr A 2 ' 1 ' 



\ m s 1 m? 

l 3/2 J 

where we have used ([5]). If ms is smaller than 2M, S decays into the gauge sector through 
the messenger loop d . The relevant interactions are 



An M 



1 i k F - 

— _ C 1 (*)M I/ _i_ _ e \*i>(xr p{i) p{i) _ _ s \( i )-p T \(0 



+ h.c, (20) 



where we neglected terms with higher orders of k(S)/M. In particular, S decays into the 
gluons and gluinos, and the decay rates are 



and 



IT] 



2 2 3 

9 ~ 64tt3M 2 { ' 



k 2 m\mQ 

T * ^ (22) 



respectively. Note that the decay rate into the gluinos is smaller than that into the gluons, 
if ms is much larger than the gluino mass [lot, 



m 5 > — m 3 , (23) 
a 3 



c Depending on the details of the UV theory, it is possible that the position of S is larger than A and the 
system cannot be described by the model given by {1} and ([2|) even for H < A. Then, the abundance of 
S will generically become larger than (fT8|) , and so, our estimate is conservative. 

d The decay into the scalars does not change our argument significantly. 
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and vice versa. Using (|T9j) . (1211) and (I22l) . we obtain the branching ratio of the gravitino 
production, 



B 



1 



3/2 



1 + r 



with 



777.3 A 



32 /4tt\ 2 
+ ~9 



7773 A 
m 3/2 M f 



(24) 



(25) 



3 \m 3 /2M P/ 

Therefore, the SUSY breaking field will dominantly decay into the gravitinos if r < 1, or 
equivalently, 

A < 2 x 10 14 GeV (*) 1 (J*L) ' . (26) 

VO.iy VlTeW UGeW v ; 



We will assume this inequality is met in the following analysis. The decay temperature of 
S is given by 



4 GeV 



r 3/2 M f 



( 9*\ 3 

uoo/ 



m 3 /2 \ : 
lGevJ 



A 



(27) 



10 14 GeV, 

where g* counts the relativistic degrees of freedom at the decay. 

Now we can estimate the gravitino abundance. Since S dominantly decays into a pair of 
the gravitinos, the gravitino abundance is given by 



3/2 



A 



Wl) vlTeVy vlGeVy vl0 8 GeVy Uo 14 GeV 



(28) 



The density parameter of the gravitino is 



0.06 k 



il a 3 \ 2 ( 7713 \ 2 / m 3 / 2 \ 



0.1 / vlTevJ vlGevJ \10 8 Gev) (l0 14 GeV 



(29) 



where h is the present Hubble parameter in units of 100 km/s/Mpc. Therefore, a right 
amount of the gravitinos can be produced by the decay of the SUSY breaking field, S. Note 
that, for the non-thermally produced gravitinos to be a dominant component of DM, k 
should not be suppressed. The reason is as follows. For a small k, the messenger mass is 
also small to keep the size of the soft masses in the SSM sector (see (181). Since the shift of 
the S field is proportional to the breaking of U(1)r symmetry, i.e., the messenger mass, the 
gravitino abundance is suppressed for a small value of k. Furthermore, depending on the 



reheating temperature and the mass spectrum of the SSM particles, the t 



nerma 



of the gravitinos and the NLSP decay may also give sizable contributions 11, 



12 



production 



13 



14 



15| 



Note also that the relatively high reheating temperature is favored, which may accommodate 



the thermal leptogenesis scenario 



Finally, let us estimate the free streaming length of the gravitinos produced by the decay 
of S. The comoving free streaming length Xps at matter-radiation equality is defined by 

Jt D a(t) 

where a(t) is the scale factor, and t D and t eq (~ 2 x 10 12 sec) denote the time at the S decay 
and at matter-radiation equality, respectively U3/2 is the velocity of the gravitino, given by 



V3/2(t) = ^ * ■ 2 UW/ (31) 
V m 3/2+ 4 \a(t)) 

where we have approximated mj ^> 07,3/2, and Op is the scale factor at the decay of S. 
Integrating (1301) yields 

\ FS ^ 7— ^=X- 1 sinh" 1 X, 

with 

= 2m 3/2 q eq 
m s a D ' 



10 6 



(lOo) (lGev) (lO^Gev) ' (33) 



where H is the Hubble parameter at present, and z eq and the red-shift and the 

scale factor at the matter-radiation equality. In the second equation of fl32|) . we have used 
Hq 1 ~ 4 x 10 3 Mpc and z eq ~ 3000. Thus, the free streaming length Xps is expressed in 
terms of the gravitino mass and the scale A, and it can be as small as 1 kpc. Interestingly, the 
recent observations on the dSph galaxies seem to exhibit a sharp cut-off around 100 pc 
in the smallest size of the galaxies, which may be explained by DM with free streaming 
length of O(100 pc) c . Our scenario may be supported by further observations in the near 
future. 



It is also possible to produce a right amount of the gravitino DM from the inflaton decay [18( , and the 
gravitino can have right free streaming length to explain the cut-off in the smallest size of the galaxies. 
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IV. CONCLUSIONS AND DISCUSSION 



The metastable SUSY breaking vacua provide a drastically simplified scheme of the gauge 
mediation. The models possess an accidental U(1)r symmetry, which is broken in the 
messenger sector. The breaking induces a linear term of the SUSY breaking field in the 
Kahler potential. We have pointed out that the SUSY breaking field is forced away from 
its minimum due to this linear term while the inflaton field dominates the energy of the 
universe. Then the gravitinos are produced when it decays, and we have shown that a right 
abundance of the gravitino DM can be realized for certain parameters. Further, the free 
streaming length of the gravitino may explain the recent observations on the smallest size 
of the dSph galaxies. 
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